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SUBJECT: BIWEEKLY REPORT, NOVEMBER 16, 1953 

To; Jay W, Forrester 

From: Scientific and Engineering Computation Group 

1. MATHEtlATICS, CODING AM) APPLICATIONS 

1»1 Introduction 

During the period covered by this report 208 coded programs were run on 
the time allocated to the Scientific and Engineering Computation (SfiEC) Group. 
These programs represent part of the work that has been carried on in 24 of the 
prolDlems that have been accepted by the S&EC Group* Progress on each of these 
problems is given below in terms of programming hours, minutes of computer time, 
and progress reports as sutmitted by the programmers in question^ 

Three new problems were initiated during this peniod* Problem #152 is a 
study of the diffusion of barium atoms in an oxide-coated cathode. This problem 
is being carried out in cooperation with B. Frost of the Digital Computer Lab- 
oratory. Problem #1^ is being carried out ty G.P. Dinneen of Lincoln Laboratory* 
This problem investigates digital methods of detecting signal from noise. The 
third problem (#155) is concerned with utilizing pressure patterns for temperature 
prediction. This work is being carried out by R. Miller of the MIT Meteorology 
Department. 

1.2 Programs and Computer Operation 

The f ollowi3[ig summary is included as a guide for interpreting the abbreviation-f^s*- 
used below. A more detailed description of the terms involved can be found in 
M-2W* 

a. The upper case letter following the problem number has the following 
significance: 

A implies the problem is NOT for academic credit, is UNsponsored. 

B implies the problem IS for academic credit, is UKsponsored. 
C implies the problem is NOT for academic credit, ^ sponsored. 
D implies the problem IS for academic credit, IS sponsored* 

The absence of a letter indicates that it is an internal SiSsEC problem. 

b. Die denotes the Division of Industrial Cooperation. 
DCL denotes the Digital Computer Laboratory* 

CMMG denotes the Cesnmittee on Machine Methods of Computation. 
DDL denotes the Division of Defense Laboratories. 
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100» Comprehensive System of Service Raatines > developed by the S&EC Group at the 
Digital Computer Laboratory for the input conversion of suitably prepared 
punched paper tapes ♦ When so requested, these routines automatically 
provide a program with suitable programmed alrithmetic, cycle-counting, 
and output facilities. ' 
:DCL Staff: Arden, 7 hours; Combelic, 6»5 hours; Demur jian, 9 hours; Denman, 
30 hours; Frankovich, 18 hours; Hazel, 7«5 hours; Helwig, 60 hours; Kopley, 
.5 hours; Porter, 20 hours; WWI,236 minutes 

A tentative version of a revised 24,6 programmed arithmetic subroutine has 
been written and is now under test as part of the comprehensive system* The new 
PA is designed primarily to correct several shortcomings of the previous PA. 
However, greater speed and increased use of mistake anticipation have been in- 
corporated into the new PA. A short written description of the subroutine will be 
issued as soon as the programs are definitely decided upon. 

Helwig 

The flad table routine for delayed print was tested in conjunction with 
J. Frankovich' s test for unassigned flads. A nxjmber of changes were made to 
improve the form, and the routine will soon be incorporated into CS. 

The program for the interpreted instructions post-mortem routine has been 
moved to occupy the end of the second bank of high-speed storage. This means that 
the programmer may now have the interpreted instructions stored in the range 40(32) 
through 3310(1736) printed out. 

Hazel 

The changes and corrections in the comprehensive system previously described 
have not yet taken effect. These revisions will take effect as soon as the necessary 
programs are recorded on unit #0. The format and content of the flad table which 
will be printed out are being changed to incorporate suggestions made by other staff 
members, 

A (24,6) arc sin subroutine has been written, tested, and included in the 
library of subroutines. 

Frankovich 

Further revisions have been made in Chapters 5 through 8 of GS Manual I. 
It is hoped to make copies of these four chapters available during the next bi- 
weekly period. 

The completion of these eight chapters will mean that we can furnish an 
applicant having no previous programming escperiemce with sufficient information 
so that he can program his problem making use of the comprehensive system. 

Denman, Kopley, Porter 

The purpose of CS Manual III is to record the inner workings of the CS 
system so that they will be comprehensible to S&EC Staff members who were not 
among the group that created the system. Annotated program manuscripts of every 
routine of the CS,(i.e., all routines recorded on MT #0) and of those control 
routines essential to its operation, will be included in the Manual. It is intend- 
ed that OS Manual III will completely document the CS system. The entire system 
is covered by the annotated program manuscripts, the text serving as a supplement. 
Work on the first draft of the text will end as of January 3I, 1953 — those 
details not covered therein being left to the manuscripts. 
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As of this date, the introduction, the First Pass Program section, and 
the Second Pass Program section have been ccmpleted. Due to current work on 
the PA, the Third Pass Program section will be ''tackled next and the PA, OT, 
Flad and Title Display, and Control programs will follow, 

Vanderburgh 

101 C, Optical Properties of Thin Metal Films on transparent backings are deter- 
mined and printed out autcjmatically by this program; the input data 
consist of the observed reflection and transmission coefficients, the 
index of the backing, the wavelength^ and the sample thickness* The 
program calculates ty means of an iterative procedure and prints out 
the index of refraction and the absorption coefficient of the film, 
the rate of variation of these constants with reflection and trans- 
mission, and the film's conductivity and dielectric constant. 
:for Professor L» Harris, Chemistry Department 
sby Dr. A> L> Loeb >(DIC^.2Q hours; J. Richmond(DIG) , 15 hours 
:DCL Staff: Denman, ^ hours; WI, m minutes 

In order to find out how fast the program converges a timed rxm. was per- 
formed on one set of data with different first estimates. The correct regult was 
n = 0*302, k = 3»269j a first estimate of n = 0.3, k » 3.3 required 10 seconds ^ 
for the film-incidence calculation with 6 additional seconds for backing incidence, 
while a first estimate n « 0.45, k = 4^5 required 18 seconds for film inctd^nce^ 
10 seconds for backing incidence. This indicates very satisfactory convergence. 

To save computer time the program has been modified for the delayed printer 
instead of the direct one. This program worked, but the column-headings program 
was accidentally omitted from the tape. 

106 G. MIT Seismic Project is concerned with the development of methods for 

locating deep reflections from underground strata in seismic pros- 
pecting. The basic method is one of prediction by means of an 
optimum linear operator, 
sfor Professor P.M. Hurley, Geology and Geophysics jProfessor G. Wadsworth, 
Mathematics Department 

I by E.A .Robins on ( Res . Assoc.), H.Briscoe, 36 hours; S. Simpson, 60 hours; 
W. Walsh, 40 hours 
:DCL,232 minutes 

Attempts to solve 9x9 and 7x7 matrices using the library subroutine 
for Craig's method were all unsuccessful. Both gave overflows while computing 
a directional vector. 

A new prediction program utilizing a graph-plotting subroutine is being 
tested but has been plagued with typing errors. 

Four frequency spectra were computed by a program using the method of 
J. Tukey. 
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108 C « An Interpretive Prograja is being developed that will accept algebraic 

equations, differential eqaations, etc* expressed on Flexowriter 
punched paper tape in ordinary n^athematical notation (within certain 
limits imposed by the Flexowriter) as input and automatically provide 
the desired solution* 
:for Dr» J*H,Laning» Jr., Instrumentation Laboratoiy 
:by J*H.Laning Jr,(pIC),25 hours; N* Zierler (DIC)« 40 hours 
■:DCL Staffs Hazel, 1 hour; WWI,' 80 minutes 

The program was successfully used for the first time during this period* 
A number of Fourier transforms were computed for problem 137 and the results are 
reported in more detail under that problem ntimber. 

Some time was lost this period because the ccmputer failed to stop ot| si 1 
on several occasions even though the stop on si 1 switch was on* 

There is an error in the part of the program that handles the printing of 
superscript variables with variable index, but it will be corrected in the next 
few days* Tapes of the program for the automatic solution of differential equations 
are now being prepared and will be tested in the coming biweekly period. 

109 C. An Airplane Pursuit-Course Program is being developed which will take account 

of airplane dynamics and projectile ballistics and thus determine an 
airplane pursuit course in three dimensions* The problem consists essent- 
ially of solving 14 simultaneous non-linear differential equations by 
the Runge-Kutta Method which is of fourth-order accuracy* 

J for Mr* J*B*Feldman, Instrumentation Laboratory 

thy M.H*Hellnian (DIGk 80 hours 

sDCL Staffs Porter, 3 hours; WWI, 24 minutes 

The two-dimensional pursuit -course problem was programmed using the GS 
before attempting the more complicated three-dimensional problem because answers 
for a particular solution of the former were available for checking purposes* 

113 C* A Stress Analysis of an L-shaped Homogenous Planar Structure is being made 
for the case of a concentrated static load* This structure is approx- 
imated by a framework of bars which will deform in the same manner as 
the prototype* This framework is then analyzed using the principles 
of virtual work and Southwell relaxation techniques. Boundary con- 
ditions have been specified for the edge of the framework so that the 
deformations of the model will conform to the actual deformations of 
the structure* 

sfor Professor J *S« Archer, Department of Civil and Sanitary Engineering 

.'by S« S ydne y (Res* Assist.CMMC) ,60 hours 

sBGL Staff J Kopley, 2 hours; WWI, 104 minutes 

A major portion of the program has been written, and some of the tapes have 
already been successfully tested* Before the program can be completely run, 
several additional tapes must be tested* 
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119 C. Spherical Wave Propagation produced by the sudden release of a spherical 
distribution of ccMpressed air in the atmosphere is being studied by 
numerical means. This involves replacing a set of non-linear hyper- 
bolic partial differential equations in 2 independent and 2 dependent 
variables by a set of difference equations written along characteristics. 
An iterative procedure is used to solve these equations* 

:for Professor C.C.Lin, Mathematics Department 

Jhy A. Ralston . CMC, 12 hours 

jDCL Staff : Gombelic, 1 hour; WWI,123 minutes 

The initial density distribution has been changed so that now at the 
center of the sphere the initial density is six times atmospheric density rather 
than the initial value of three tijnes atmospheric density used in the previous work. 

Two programs for different initial mesh widths were run using the new initial 
density distribution. The results from both indicate the formation of a shock this 
time. With the formation of the shock various difficulties arise in connection with 
the numerical process being used. These difficulties necessitate a great deal of 
work off the computer before the computation can^ continued on WWI. 

123 C. Earth Resistivity measurements are used to calculate the Slichter kernel 
function which, in special cases, can be analyzed to give the actual 
distribution of resistivity. The method involves least-square fitting 
a set of polynomials to the measured surface-potential function and 
integrating the product of this set and the zero-order Bessel function* 

:for P.M.Hurley, Department of Geology and Geophysics, DIG 5-6915 

*^ K. Vozoff (Res> Assist. CMMC),4 hours 

:DGL: WWI, 34 minutes 

The tape for computing Jq(x) and J-j_(x) was completed. The program, written 
in CS, computes the two functions simultaneously using a difference equation. This 
corrected subprogram will be inserted in T-2698, the program for computing the 
Slichter kernel, and the latter will be retested. 

126 G. A Data Reduction Program for use in the Servomechanisms Laboratory is 
being developed in separate stages to be combined at a later date. 
The first stage is concerned with devising a program to fit poly^ 
nomials to arbitrary empirical functions using a <least squared error 
criterion. The procedure makes use of Legendre polynomials and matrix 
multiplication. 
:for J.E.Ward, Servomechanisms Laboratory, DIG No. 7138, AF33( 616)2038 
:by D.T.Ross (DIG).60 hours; D. Hamilton, 60 hours; R. Turyn(DIG),80 hours 
:DCLi WWI, 196 minutes ^r 

The Mistake Diagnosis Routine (MDR) program has now been put in final form, 
and several applications have been made to the Polynomial Fit Program for the 
investigation of round-off errors. Copying errors in the final version of the 
MDR and in a small routine necessary to sample one section of the Polynomial Fit 
Program have delayed results. A check-order alarm during the execution of a bo 
order on the latest run is unexplained. The report on the MDR is being published. 

The Lagrange Interpolation and P^st-Mortem programs have not yet operated 
successfully. Modifications are being made and testing continues. 

Excellent results were obtained by fitting a fifth-order polynomial to an 
empirical function using the Polynomial Fit. The maximum error was 0.16 out of 
lil',500 with variance of all the errors •0015. 
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131 • The Training of New Personnels Tours and Demonstrations are among those 

activities included in this problem* Generally speaking, any approved 
staff problem relating to training and/or demonstrations is considered 
to be in this categoiy. 
jDCL Staff: Koplev . 3.5 hours; WWI, 52 minutes 

Three tours of WWI, listed below, were given during the past biweekly period. 
All three demonstrated several programs on the computer and included a general tour 
of the installation. The first two groups mentioned below were also given a Flexo- 
writer demonstration. 

Date Time Nuinber Represented 



November 3 5-6 PM 23 19 from Sylvania 

4 self 

November Ik 10-11 AH 10 Institute for Radio Engineers 

November 15 8-9 FM 13 Old South Ghiirch, Reading, Mass. 

132 C» Subroutines for the Numerically Controlled Milling; Machine are being 

revised, and tested, Therset of subroutines facilitates programming 
of the computations involved in the preparation of numerical data used 
to control the milling machine. The subroutines involve routine num- 
erical and logical operations* 

sfor J . 0,McDonough, Servomechanisms Laboratory, DIG Project 6873 

•^y J«H,Runvon (E.E^Res^ Assist,), 10 hours 

;DCL Staffs Combeiic, 2 hours; "WWI * 8 minutes 

Testing of subroutines for finding points on the series I6 symmetrical 
section and camber line was continued. One error, an incorrect address, was found 
and corrected in the routine for the symmetrical section* The routine then worked, 
properly. The camber line routine still has not operated properly* 

A 90 section of a supersonic nozzle was cut using the tape prepared on 
WWI, Results were very satisfactory, 

137 D* Investigation of Atmospheric Turbulence as a noise input to airborne 
control systems, A stationary random process is assumed so that 
the methods of generalized harmonic analysis may be used to describe 
the turbulence components in terms of their power spectral densities* 
:for Professor R,C*Seamans .Department of Aeronautical Engineering 
iW R, A, Summers (Res, Assist,). 8 hours; N.Zierler,^ hours ;C .Block, 2 hours 
sDCL: 'WWI, 87 minutes 

The generalized program of problem 108 has been successfully utilized by 
N, Zierler to compute a Fourier transform( power spectral density) , The results 
show all positive "power" out to 10 radians /second, except for one isolated and 
very small negative number at 4 radians /second. This result is also in agreement 
with an independent GPC computation. 

It must therefore be concluded that the Fourier Transform program of 
problem 107, which (for the same data as above) produced all negative power above 
5 rad,/sec,, contains some inherent innaccuracy, perhaps in the approximate recursion 
formulae for sine and cosine. This program is being checked with a v^ew to 
eventually using it, in its corrected form ,for the Fourier transform calculations | 
the generalized program of problem 108, while being very la^ccurate and reliable, may 
not be fast enough for taking the large number of F.ourier transforms required in 



M-2524 page 7 

problem 137 • 

Completion of problem 137 depends on a satisfactory correction of the 
Fourier transform program of problem 107. 

138 B. Spheroidal Wave Functions are solutions of the scalar Helmholtz equation 
separated in spheroidal coordinates. A program has been developed 
for tabulating both the coefficients obtained by expanding the angular 
solutions of the first kind in associated Legendre functions and the 
coefficients obtained "by expanding the radial solutions of the first 
kind in spherical Bessel functions. By analytically substituting these 
expansions in the separated ordinary differential equations, 3-'t'erm 
recursion relations are obtained for the coef f icients^, and the radial" 
and angular coefficients are found to be simply related. Both sets of 
coefficients are then determined by applying an iterative procedure to 
a continued fraction equation derived from one of the 3-"term recursion 
relations. The iteration proceeds until a value of the \inknown separ- 
ation constant of the differential equation which appears in the con- 
tinued fraction equation is found such that it makes the coefficients 
compatible with their boundary conditions. 

:for Professor P.M.Morse, Physics Department 

;by F.J. Corbato (Res.Assist.CMMG) ,30 hours: J.D.C.Little CRes.Assist. CMMC), 
30 hours 

:DCL Staff: Combelic, 2.5 hours; WWI* 164 minutes 

Production was barely started when an error was discovered in the present 
PA buffer. Statistically the buffer contents are incorrect one in every 1024 its 
operations, "'^his occurs because a round- off made in the 25th binary digit occasion- 
ally causes a change in the first 15 "significant binary digits. ' The 
effect of this error is to reduce the accuracy from 30 to 15 binary digits. J-his 
mistake will be eliminated in the new CS PA being written by Heiiwj.gi. 

1^* Summer Session System consists of a conversion program, an interpretive 
routine, and mistake diagnostic routines stored in WWI. A special 
mnemonic instruction code has been developed for use with this system 
thus simulating a computer with characteristics quite different from 
those of WWI. This Sujumer Session (SS) ccmputer was developed for the 
use of students participating in the MIT 1953 summer session course on 
"Digital Computers and Their Applications" . The SS computer is being 
used in the E,E. Department course 6*535 and is available to programmers 
with suitable problems. 
:DCL Staff: Combe3.ic,~ 42.5:hour5J Fr^nkovich, 4;hours; Helyig, 10 hours 1 
Hoy, 58 hours; Siegel . 45 hours ;^WWI, 184 minutes 

During this biweekly period, attention was principally directed toward 
developing the automatic title display described in the preceding biweekly report. 
A program for the title display on the oscilloscope, which is the most difficult 
case of the three, is undergoing test. Satisfactory operation has not yet been 
attained, but it is expected that the scope-title feature will be in operation 
before the end of the next biweekly period. 

Several programs were run on the SS computer for students in Course 6«535« 
The operation of the computer was entirely satisfactory. 

A significant mistake in the conversion of the floating-point number zero 
was located and corrected. 
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141* S&EG Subroutine Study has been undertaken for the final testing of sub- 
routines selected for incorporation into the Library of Subroa tines. 
Although very little effort is going into the specific development of 
subroutines , programs that have been written for other S&EC problems 
and seem to be of general use will be suitably modified for the S&EC 
Library ♦ 
:DCL Staff: Arden . 6 hours; Dennian,10 hours; WWI, 9 minutes 

A routine making simple formats available has been added to 2756m5 and is 
currently available in the Library of Subroutines as 3217m5. 

1^3 '3* • The Vibrational Frequency Spectrijm of a Copper Crystal is to be determined 
by solving a 3 x 3 secular determinant, each terra of which consists of 
a finite Fourier Series of 12 terms. This equation must be solved for 
24, 495 different ^values of the wave-propagation vector,, 

:for Professors B^E.Warren and J,C^Slater, Physics Department 

*^y E>HsJacobsen (Res. Assist*), 12 hours- 

: DCL Staff; Conibelic, 5 hours; WWI, 205 minutes 

A double integration concerned with describing the scattering of X-rays 
from a crystal lattice has been performed successfully 21 times (i.e., for 21 
different values of a parameter). The double integration must be performed 6 
more times to complete this phase of the problem. 

144 C. 3,elf-Consistent Molecular Qrbitals are the optimum choices of linear 

combinations of atomic orbitals determined through a process described 
as a self -consistent field approximation. The numerical procedure involve 
matrix -vector multiplications, vector additions, and matrix diagonal- 
ization, 

J for Professor J „G.» Slater, Physics Department 

J by Dr. A. Meckler (DDL). 5 hours 

sDCL Staffs Arden, 8 hours; WWI, 83 minutes 

The program has run through one iteration without an alarm. The next 
problem to be investigated is the convergence behavior of the entire scheme. 
With different initial guesses to the density matrix, the machine will be forced 
through different iteration routes, revealing the possibility of divergence and 
the rate of convergence, 

14? C, Energy Bands in Crystals are being studied hy finding solutions of the 
corresponding second order linear differential equation satisfying 
"boundary conditions at the origin. The solutions are found approx- 
imately by using the Gauss-Jackson formula for forward integration. 
The solutions and their first derivatives are to be combined in a 
sum, the weighting factors being functions of an independent para- 
meter. 

J for Professor J. C. Slater, Physics Department,, DIG No. 6853 

iby Dr. D.J.Howarth .(DIC). 16 hours 

:DCL Staff J Arden, 12.5 hours; Porter, 2 hours; WWI, 109 minutes 

The recurrence relation which was originally used to ccanpute the spherical 
Bessel functions j*(x) = J^ Jj ,i(x) proved to give a large accumulation of error 

2l -1 ; "^ 

when -~ > li It has therefore been decided to compute j« (x) as follows: 

defining Av (x) ^P(v+l) 2^^ x"'^J^ (x),Aj,(x) varies slowly with x, andA^(0)=l, 
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For high y»Av ^^y ^® obtained from a power series, and use can then be made of 
the backward recurrence relation 

This proceeds without loss of accuracy provided /. Y\j+i^ ^ ^* There is insig- 
nificant loss of accuracy for v^f for all x used here (x<8), A subroutine 
to perform the above calculation has been programmed. The testing of the final 
part of the main program has been completed; some difficult errors have been 
eliminated. 

149 C . Digital Methods of Detecting Signal from Moise are being investigated. 

A sequence of binary numbers will simulate the message wherein regions 
of high density of ones are signal regions and those with low density 
of ones are noise regions. Various methods of detecting the change 
from one region, to another, as well as the length and midpoint of the 
signal regions are being studied. 

:for J.V.Harrington, Lincoln Laboratory 

J by G.P.Dinneen, Lincoln Laboratory 

I DCLi^ ¥WI, 14 minutes 

Two basic schemes are to be investigated Jfirst. One is a sequential 
observer as described in Lincoln Laboratory Technical Report No. 20 and the other 
a success run observer. The Sequential Observer employs a biased counter with an 
upper and lower threshold and another counter which measures target length, 
triggered by the attainment of the upper threshold and reset when the lower thres- 
hold is reached* The success run has a single counter which counts every pulse 
and is reset by some configuration of pulses. These two schemes or modifications 
of them will be tried for various signal to noise ratios and signal lengths* 

A program for two sequential observers and a* success -run observer was run 
off, and results obtained. The success-^run observer appeared more selective for 
this case. However, it was determined that the number of bits per experiment 
should be increased and a new method of handling the end of the experiment should 
be devised. These changes were incorporated in a modified program whi(5h should 
be run during the next period. 

152 B . Diffusion in an Oxide-Coated Cathode is a program to calculate the effects 
of combined thermal and electrolytic diffusion that occTor in an oxide- 
coated cathode when current is caused to flow through the cathode. 
2 for W. B, Nottingham, Pliysics Department, DIG No. 63^5 
iby H.B.rrost( Res.Assist. E.E.Dept.), 8 hours 
iDGL Staff; Denman, 2 hours j Porter, 6 hours; WWI, 336 minutes 

The equation governing the process is a second order non-linear partial 
differential equation of parabolic form, as follows: 

where n is the concentration of bari\am donor centers, D is a thermal diffusion 
constant, and C is a parameter depending upon coating thickness, the coating 
conductivity, the current through the cathode, and temperature. The concentration 
n is normalized and has the value 1 for G<x<a, ts=0. The basic boundary condit- 
ion, resulting from the fact that donor centers do not migrate across the 
boTondary. is 
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an aaixiliaiy condition, very useful as a check on errors in computation is 
— r^ ( ^ n(x,t)dx ) = 



o 

This .problem may be Bet up in the following difference equation form: 



^(x,t4-l) -^.t ^ ^^^, sl Vl^t ">PSTt ^ Vl.t "^^.tTVl.t ) 



\ ^x h^ 

X 

for 0<x<N* 

2 
Setting h.«^ —x , Y= V and h = f 
t 2D fl "Y^ ^ 

yields 






The boundary oonditions may be approximated to the second order at constant t by 

if! 

and "~~~ 



%-^n.L^.^(ii)'-im"ri".-(i>'y 



There are two results of prime interest which have physical significance. These are 
^ ^ J \[n(x,t) 



and f^Ct) = ^n(0,t) 

The integral may be approximated satsifactorily by trapezoidal integration. 

During the past period machine runs have been made to eliminate programming 
difficulties. With a good program, results have been obtained for ^ *== 0«03 using 
approximately 300 minutes of machine time. The restilts can be compared to steady— 
state solutions available by analytic methods. The comparison shows a truncation 
error in the integral 

n(x) dx 

of O.OOIS/^, with a probability that this error would increase further if the 
problem were run longer. 
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Efforts are being directed toward reducing the running time and improving 
the accuracy. Most of the machine time is used for calculating square roots. 
A considerabl,e saving can be accomplished through the use of a square-root sub- 
routine which makes use of available good approximations to the roots desired. 
Such a program has been written, 

A better approximation to the boundary rconditions has been derived using 
third differences. This approximation has been tried for the first three steps 
of calculation, with a reduction in the total truncation error by greater than 
a factor of two and a reduction in the rate of increase at the third point hp a 
factor of five. This approximation looks ^ery favorable. It requires less 
computation than the present boundary conditions since it requires one less square 
root operation, 

155 B« Synoptic Climatology^ A multiple regression formula is used to predict 
temperatures from pressure distributions described by Tschebycheff 
polynomials* The matrix of scalar products which is used in the 
calculation of the coefficients of the multiple-regression system is 
being calculated on "WWI, 

sfor Professor T.F.Malone, Meteorology Department 

:by R. Miller (DIG) 

sDCL Staff t Arden, 10 hours; Wl, 26 minutes 

Several errors in this program have been eliminated* 



1.3 Operating Statistics 

Computer Time 

The following indicates the distribution of WWI time allocated to the 
S&EC Group* 



Programs 
Conversion 
Magnetic Di^um Test 
Magnetic Tape Test 
Scope Calibration 
Demonstrations (#131) 
Total Time Used 
Total Time Assigned 
Usable Time, Percentage 
Number of Programs 



38 hours, 38 
17 hours, 56 
20 
16 
78 
H6 



minutes 
minutes 
minutes 
minutes 
minutes 
minutes 



59 hours, 14 

60 hours, 55 
97.67i 

208 



minutes 
minutes 



1*4 Summary of Tape Room Bulletin Board Memoranda (I. Hazel) 

(These memos are intended to inform programmers of changes in coding 
procedure, Wfl facilities, etc.) 

The Cycle Count Instructions 

For a detailed description of the Cycle Count Instructions, see Bulletin 
Board Memo #44, which is posted in the Tape Room, Barta 110. 
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Multiple-Buffer Subroutines 

This subroutine actually modifies the PA routine so that the PA 
routine refers to 2 constants in the multiple-buffer subroutine ♦ The reference 
to these two registers is made during the executi'Dn of each interpretive . 
instruction whether or not the instruction refers to the multiple buffers • For 
this reason, the multiple -buffer subroutine - must remain pennanently in high- 
speed storage at a permanent location for a given program. 

The use of the regular buffer of the PA (symbolized by b) will have no 
meaning when the multiple -buffer subroutine is used. 



M-2524 pg^gg 2,J 

2. GCHPUTER ENGINEERING 

2^1 WWI System Operation 

2«11 Core Memory (L.L.Holmes, R.J.Roberts) 

The majority of low margins in the system have been cleared up. Excursions 
have been set in on the marginal-checking lines for core-storage control, and all 
margins are reasonable. 

The installation of the parity-check system for the drum was delayed until 
November 16. 

A program for generating a "worst possible pattern" in the cores should be 
ready for trial in a few days. 

2.12 Atixiliary Drum (K.E.McVicar) 

With the exception of one instance the auxiliary drum has operated without 
known error during the last biweekly period. On one occasion an intermittent 
readout was obtained during program operation. The trouble was traced to writing 
between the slots on the dram surface. A visual check had been made of the dnon 
tracks on the previous day, and no serious defects were noted in the waveforms. 
It is assTimed that the writing between slots was caused by the marginal checking 
which was done on the morning of the day the readout failure occurred. Investi- 
gation of the problem revealed that three lines could cause writing between the 
slots in case of certain failures caused hy varying the operating voltages. 
These lines have been temporarily removed from the routine marginal-checking pro- 
cedures* In the meantime work is in progress on a new program for marginal check- 
ing the drum on the read operation only, except where the write function itself 
is to be checked* 

2.13 Typewriter and Paper Tape (L.H.Norcott) 

During the past two weeks we cempleted routine overhaul on the remainder of 
our long-carriage Flexowriters. . 

I have been working with Farnsworth and Perry to locate the cause of inter- 
mittent errors appearing recently in both tapes and printed copy produced on the 
delayed- output Flexowriters. The work is not completed, but it now appears that 
the trouble is in the magnetic-tape equipment rather than in the Flexowriters* 

2.2 Terminal Equipment 

2e21 Magnetic Tape (E .P. Farnsworth) 

Several defective tubes, crystals, and tapes were located by marginal check- 
ing and replaced during this period. Other troubles experienced were caused by 
a defective paper-tape punch on the delayed-output Flexowriter, and failure of 
programmers to use the si stop magnetic tape order when recording. 

Minnesota Mining is now working on a triacetate tape for us to replace mylar 
as they have been unable to get a definite commitment frcm Dupont on the supply 
of unspliced mylar. Triacetate is stronger than acetate, though not as strong as 
'iTQrlar, but it does not stretch like mylar. Stretching has been the main cause of 
failure of mylar tape. 
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3. LIBRARY ACCESSIONS LIST 

The following material has been received in the Library, ¥2-325 • 
Library Files 
No, Source ^ . Title_ 



2566 Willow Run Rsch^ Ctr« 
257G Oim/London 

2577 Lincoln Laboratory 



MIDAC Automatic Computer (Brochure) 

Computers at the Federal Institute of 
Technology, Zurich 

Symbolic Design of Digital Ccanputers 



Laboratory Files 

lot Author 



M-.2461 CFleming-M.Demurjian 



Title 



Instructions for the Use of the Flexowriter- 
Verifier Equipment, 



M-2465 



S&EC Biweekly, 10-19-53 
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^^ AIMINISTRATION AM) PERSOMEL 

Terminated Staff ( J .C .Proctor) 

G, A# Murdoch 
D* Goldenberg 

New NoH'-Staff (R. A. Osborne) 

James Mahoney has retmmed to the Laboratory and is now a member of the 
T/?hittemore BmUding janitor crew» 

Charles Muhle is a new technician in Group 64, temporarily assigned to the 
Constniction Shop, 

John A. G*Brien is a new technician in the Construction Shop. 

Ralph Porter is another new technician in the Construction Shop. 

Terminated Hon^^taff 

Gordon Kane 
Michele Rhetme 
James MacDonald 



